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We connect between the problem of thermodynamics of localized magnetic moments in a Dirac 
semimetal, the interaction with relativistic electrons leading to the effective ferromagnetic exchange 
between the moments, and the existing theories dealing with long-range exchange interaction. We 
point out that the results of high-temperature expansion for the free energy of a dilute ensemble of 
magnetic impurities in the semimetal performed by V. Cheianov et al. (Phys. Rev. B 86, 054424 
(2012)) give an indication to the existence of a new disordered fixed point in such model. 
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Long-distance exchange interaction between mag¬ 
netic moments mediated by the mobile carriers is 
known as the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
exchange^. Recently there appeared interest in RKKY 
interaction in the class of materials in which the low- 
energy electron excitations resemble massless Dirac par¬ 
ticles: graphene^i^, chiral metals formed at the surface of 
topological insulators^i^, and silicene^. There is a pecu¬ 
liarity of the RKKY exchange in such conductors which 
make them qualitatively different from usual metals: the 
Friedel oscillations are either absent or commensurate 
with the lattice^. In particular, quite a few papers stud¬ 
ied collective behavior of magnetic adatoms randomly 
distributed on the surface of a topological insulator— 

Our communication is inspired by the very interesting 
publication by Cheianov et al— The authors considered 
the high-temperature expansion in the inverse tempera¬ 
ture for the disorder averaged magnetic susceptibility of a 
dilute ensemble of Ising magnetic impurities in a 2d Dirac 
semimetal. From this expansion they found the critical 
temperature of the ferromagnetic phase transition and 
the magnetic susceptibility critical exponent. 

We suggest that high-temperature expansions can 
shed light on critical behavior of the long-range exchange 
interaction models, both with and without quenched dis¬ 
order. That is why we decided to briefly sum up vast 
body of existing results in the field, obtained by the 
renormalization group (RG) analysis, numerical simula¬ 
tions, etc, emphasising still open problems in the theory. 

RKKY effective exchange interaction between a pair of 
localized magnetic moments described by spins Si and S 2 
has a very simple structure 

Hrkky = (1) 

where J{R) = /^x(i?), I is the exchange interaction be¬ 
tween the localized magnetic moment and itinerant elec¬ 
trons, i? is the distance between the magnetic moments, 
and 

x(i?) = -i^'^^e(R;r)e(-R;-T)dr (2) 


is the free electrons static real space spin susceptibility^^. 
The Matsubara Green’s function Q isii 

e(R,T) =-(r,c(R,T)ct(O,0)). (3) 

Further on we assume that T = 0 and the Fermi energy 
is at the Dirac points. Then the Green’s function is 


Q{R-,t) 


—sign(T)0 


ik-R-i;fc|r| 

{2'k)'^ 


(4) 


where d is the dimensionality of the space, Ul is the vol¬ 
ume of the elementary cell, and v is the velocity of elec¬ 
trons. Performing integration in Eq. we get 


G{R;t) 


e(R;r) 


sign(r)Ou|r| 

(R2 + i;2^2)3/2’ 

sign(T)DRu|r| ^ ^ ^ 

( i ?2 + i ; 2 ^ 2)2 ’ 


(5) 

( 6 ) 


Next performing integration in Eq. m we obtain 


x{R) ~ 


vR^ ’ 


d = 2 
d = 3. 


(7) 

( 8 ) 


(Actually Eqs. ([7]) and ([5]) can be obtained just from 
dimensionality considerations; numerical coefficients are 
anyhow of no interest to us.) 

The problem of thermodynamics of magnetic mo¬ 
ments, forming a periodic lattice, with an isotropic 
n—component order parameter and algebraically decay¬ 
ing ferromagnetic exchange interactions J{R) 

J(i?) - 1/7?'^+°', (9) 


corresponds to the effective 0{n) Hamiltoniani^ 



( 10 ) 
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(Transition from the discreet exchange Hamiltonian to 
Landau-Ginzburg one is explained, for example, in 
Ref.—.) For the short-range exchange 6 = 0. For the 

long-range exchange 6, c 0; the term is gener¬ 

ated dynamically even if the initial microscopic Hamilto¬ 
nian is purely long-range. 

Study of such a Hamiltonian has a very long history. 
We refer the reader to Ref— for the list of works pub¬ 
lished before 1997. The phase diagram for this model 
was proposed in an early seminal contribution^^. For 
d > min(2CT, 4) the model is characterised by Gaussian 
fixed point. In this regime for all n we have 


VG 

= 2-a 

(11) 

IG 

= 1 

(12) 

VG 

= 1/cr. 

(13) 


Eqs. (fTTll and (IT^ differ from those of Landau-Ginzburg 
theory: rj = Q, v = 1/2. This is due to the fact that for 
the Hamiltonian (1101) field correlation function in param¬ 
agnetic phase (in momentum representation and for small 
momentum) in the mean field approximation is 


G{q) 


1 

+ r ’ 


in distinction to the traditional one 


G{q) 


1 

cq'^ + r 


(14) 


(15) 


(In both cases r ^ t = {T — Tc)/Tc.) The critical expo¬ 
nent r] is definedi^ by Equation (for T = Tc) 


G{q) ^ q-^+^, (16) 


which explains Eq. m- 

For d < min(2cr, 4) two non Gaussian fixed points com¬ 
pete with each other: the Wilson-Fisher or short-range 
(SR) fixed point^S^ meaning that the model is equivalent 
to one with short-range interactions and the Fisher-Ma- 
Nickel or long-range (LR) fixed poini— , specific for long- 
range interaction. The case © lies on the boundary of 
the ’’classical” region and the long-range fixed point re¬ 
gion, where the observables differ from those of the mean 
field theory by logarithmic factors. In the critical region 
above the critical temperature Tc the correlation length, 
susceptibility and heat capacity vary 


e(r)^t-i(int-i)”' 

x(r)^t-i(lnt-i)”', (17) 

C~(lnt“^)^"' (n < 4); C ~ lnlnt“^(n = 4). 

where n' = (n-|-2)/(n+8). Thus for Ising model n' = 1/3, 
for XY model n' = 2/5, and for isotropic Heisenberg 
model n' = 5/11. For n > 4 the specific heat remains 
finite and has no jump. In the critical region below the 


critical temperature Tc the spontaneous magnetization 
varies as 


m(T) ~ (18) 


It is worth mentioning that the predicted logarithmic 
corrections were accurately observed in extensive Monte 
Garlo simulations of Ising modelsii. It would be in¬ 
teresting to see to what extent do these results corre¬ 
spond to one obtained from high-temperature expan¬ 
sions?. We wonder, whether a kind of special treatment 
of hig-temperature expansions proposed in Ref— to ex¬ 
tract logarithmic corrections can be of some help. 

Here probably a simple explanation, why the upper 
critical dimension in the model with short-range ex¬ 
change is 4, and in the model with long-range exchange 
m is 2a (for a < 2), would be relevant. Scaling transfor¬ 
mation of the Hamiltonian cni) starts with writing down 
Hamiltonian m in momentum representation (integra¬ 
tion with respect to q is limited by some ultraviolet cutoff 
A). We perform integration with respect to q, satisfying 
A/s < q < A, where s ^ 1—. In second order of pertur¬ 
bation theory in go the only graph which is necessary to 
take into account is proportional to 


/ A/s<q<A 


d‘^q 

(27r)« 


G^q), 


(19) 


where Green’s function is calculated for a = 0. Integral 
(fmi) with Green’s function (fT51) is 


/ A/s<q<A 


d’^q 

(27r)° 


G\q) 



A-^{s^ - 1) 

e 


( 20 ) 


where e = 4 — d. On the other hand, for Green’s fuc- 

tiondmi 


I A/s<q<A 


d‘^q 

{2t:)° 


G^q) 



A-'^js'^ - 1 ) 
€ 


( 21 ) 


where e = 2a — d. In both cases logarithmic dependence 
upon s 


! A/s<q<A 


d<^q 

{2nY 


G^{q) ~ Ins 


( 22 ) 


corresponds to e = 0. 

In the lowest order approximation (with respect to g), 
of two terms bq'^ and cq^ the term with the lower degree is 
relevant and the term with the higher degree is irrelevant. 
Thus the transition between the fixed points for d = 4 
corresponds to ct = 2. Fisher et al— assumed that this 
remains true for any d, and, while in the long-range fixed 
point region the exponent 7 is a nontrivial function of a 
and d, simple Eq. m is valid there for the exponent g. 

The last statement can be justified, in particular, in the 
large-n limit of 0{n) model^i, as it was done in Ref—. 
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We put in Hamiltonian (11011 c = 0, 6 = 1. and introduce 
the auxiliary imaginary field A(a;) conjugate to (fP". 


exp 


DX exp 


/ d‘^x 


/ 

2^ 8n \ J 


d^a 


^X-- 


2g 


nrX 

9 


■ (23) 


We keep the longitudinal component </>! (fixed by a van¬ 
ishing external field) and integrate on n — 1 transverse 
components of (f). We finally set 4>i = and arrive to 
a reduced Hamiltonian: 


H{X,(p) 




n 
n — 1 


Trln [-A^ 


A], 


25 . 

(24) 


The large n limit is thus given by the saddle point equa¬ 
tions in the two fields ip and A and the corrections are 
the usual loop expansion. In the absence of space varying 
external field, we obtain the equations 


ipX 


X — r — 


9 2 
2 ^ 


0 

9 f d'^q 1 

27 {2(j))W^ + y 


(25) 


The previous equations can be solved easily. At and 
above Tc the magnetization ip vanishes, at Tc A = 0, 
above X ^ 0. Thus above the saddle-point equation 
reads 


t _ 2 1 r d’^q 

A 9^ {2-kY J q<^ {q'^ + A) ’ 


(26) 


where r — Vc = ft, and t is proportional to T — T^. For 
d > 2a the integral converges when A vanishes and one 
obtains the mean field result 


^ = A-i/-"(27) 


by taking into account higher order terms in the RG cal¬ 
culations, predicted that the change of behavior from the 
intermediate to the SR regime takes place at cr = 2 — risR. 

Many other studies also have considered this problem 
of tr = 2 with various conclusions. In particular, van 
Enter— obtained that for n > 2, for the classical and 
quantum XY models long-range perturbations are rele¬ 
vant in the regime cr = 2 in contradiction with Sak re¬ 
sults. The same statement for arbitrary n was made in 
Ref—. The Sak scenario was also challenged in Ref—, 
which presents results of a Monte Carlo study for the fer¬ 
romagnetic Ising model with long-range interactions in 
two dimensions. The author claims in addition that the 
results close to the change of regime from intermediate to 
SR (cr > 2) do not agree with the renormalization group 
predictions. 

A first numerical study of the exponent iq for d = 2 as 
a function of a has already been done in Ref—. In par¬ 
ticular, the authors obtained in the intermediate regime 
{d/2 < a < 2) a result well described by the exponent 
5 = 5 g = 2 — crupto 2 — cr = rjsR and rj = rjsR for larger 
cr. In the subsequent paper—, the authors claim that 
the boundary between the SR and the LR fixed points 
for d = 2 corresponds to cr = 7/4. In a field-theoretic 
approaches it was proved, to all orders in perturbation 
theory, the stability of the SR fixed point for cr > 2 — rjsR 
and of its LR counterpart for a < 2 — r]RR, where tirr 
is the anomalous dimension of the field, evaluated at the 
long-range fixed point—. Quite recent numerical analy¬ 
sis of the problem was presented in Ref—. By including 
the subdominant power law, the numerical data are con¬ 
sistent with the standard renormalization group (RG) 
prediction by Ref—. 

These debates have a practical importance for the case 
dHl), corresponding to cr = 2. If this case is described by 
the SR fixed point, we have^S Eq. (EH) and (in quadratic 
expansion with respect to e) 

1 _ 1 _ 7 ^ + 2 \ £ _ {n + 2){7n + 20) /£ 

'ysR \n + 8)a (n-I- 8 )^ W 



i.e. V = 1/a. For a < d < 2a the integral in the r.h.s. of 
( 1 ^ diverges near as i.e. 

^ = l/(c;-cr), (28) 

and from the scaling law v{d—2+r]) = 2(3 and the relation 
(3=1/2 valid in given approximation one recovers Eq. 

m- 

Although the general outline of the phase diagraroi^ 
has been widely accepted, the location of the boundary 
between the SR and the LR fixed point has become the 
scene of a debate. The reason of the objections to the ini¬ 
tial position of such boundary at cr = 2 is simple. From 
the conjecture follows that \iina ^2 5 = 0. Together with 
this for cr > 2, d < 4 the critical exponents assume their 
SR values, with positive value of r]sR. Then it would 
imply a jump of the exponent rj from 0 up to rjsn at 
cr = 2. This contradiction was removed by Sabi^, who, 


with e = 2a — d = 1 , and 


Q{a) = a 1 /( 1 ) - 2i/ -fi V'(o-) 


(30) 


where ip is the logarithmic derivative of the gamma func¬ 
tion. If this case is described by the LR fixed point, 
we have the e-expansion for the critical exponents of the 
0{n) model,— (in the same approximation): 


ILR 

ILR 


= 1 


2(n-fi8)2 
n -I- 2 


n -I- 2 
4(n-t8)' 


■{rP + 22n + 52)e^ 


2{n + 8)' 




(31) 


with e = 4 — d = 1. High temperature expansions may 
supply important argument in the debates. 

The presence of quenched disorder can qualitatively 
change critical behavior of a magnetic system. A general 
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argument^ shows that one should expect a new type of 
critical behavior for the random system, distinct from 
that of a pure one, whenever the specific heat of the pure 
system diverges at the transition temperature. This cer¬ 
tainly happens for both cases © and ([5]). 

Fixed point 0(e) where e = 4 — d for short-range in¬ 
teraction model with weak quenched disorder for n > 1 
was found by Harris and Lubensky^^. Critical exponents 
calculated for this critical point are 


n(5n — 8) 2 

~ 256(n-l)2^ 

(32) 

3n 

(33) 

ISR - 1 17^- 

16(n — 1 ) 


For n = 1 Khmelnitskii found another disordered fixed 
point of order Critical exponents calculated for 

this critical point^i^^ are (e > 0 ) 

VSR = + (34) 

7S. = l + +0(e). (35) 

Results for magnetic susceptibility and heat capacity in 
the critical region at d = 4 can be summed up as^l 


X(T) ^ 

t ^ exp 


{int-r 

(36) 

C{T) r. 

exp 

-2 {D 


(37) 


where D = 6/53, and a = 1 / 2,7 = 0^—, a = 1.24,7 = 
_ / jilAl 

We are aware of a single paper where a model con¬ 
taining both disorder and long-range interaction was 
studied^. There critical properties of a random Ising 
model with long-range isotropic interactions m were 
analysed by using renormalisation group methods in an 
expansion in e = 2tT — d. For e > 0 the critical behaviour 
was described by a stable fixed point 0(i/e). Like in 
previous papers considering no-quenched disorder case, 
it was found that when cr = 2 — ryg/j the system crosses 
over smoothly to SR behaviour. The peculiarity is that 
for the random fixed point r/sR < 0 (see Eq. (121 )• Thus 
the crossover to SR behaviour is analysed and takes place 
when (7 = 2-1- e/106, 4 — d > 0. Thus, according to this 
work, the case (d with added weak quenched disorder 
corresponds to the LR fixed point. 

It is unclear whether fixed point obtained for weak 
disorder potential scattering added to Hamiltonian (fTUl) 
is the only possible one for models with disorder. 
We believe that critical exponent obtained by high- 
temperature expansions ii>i^ (7 = 1.4) is an indication 
of the existence of a new fixed point existing for the 
strong gas disorder model considered there. Hence ad¬ 
ditional studies of high temperature expansions in that 
model (with different dimensions, different exchange de¬ 
cay laws etc.) would be of much interest. 

The authors are grateful to D. E. Khmelnitskii for very 
illuminating discussions. 
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